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Abstract — We derive the capacity region of the Gaussian 
version of Willems's two-user MAC with conferencing encoders. 
This setting differs from the classical MAC in that, prior to each 
transmission block, the two transmitters can communicate with 
each other over noise-free bit-pipes of given capacities. 

The derivation requires a new technique for proving the 
optimality of Gaussian input distributions in certain mutual 
information maximizations under a Markov constraint. 

We also consider a Costa-type extension of the Gaussian MAC 
with conferencing encoders. In this extension, the channel can be 
described as a two-user MAC with Gaussian noise and Gaussian 
interference where the interference is known non-causally to the 
encoders but not to the decoder. We show that as in Costa's 
setting the interference sequence can be perfectly canceled, i.e., 
that the capacity region without interference can be achieved. 

I. Introduction 

We consider a communication scenario known as the MAC 
with conferencing encoders where two transmitters wish to 
transmit independent messages to a single receiver. Prior to 
each transmission block, the two encoders are allowed to 
hold a conference, i.e., they can communicate with each other 
over noise-free bit-pipes of given capacities. Special cases 
are the classical multiple-access setting, where the encoders 
are ignorant of each others messages (both bit-pipes of zero 
capacities); the fully-cooperative setting (both bit-pipes of 
infinite capacities); and the asymmetric message sets setting, 
where one of the encoders is fully cognizant of the message 
the other encoder intends to send (the pipe from the cognizant 
transmitter to the non-cognizant transmitter of zero capacity 
and the other pipe of infinite capacity). 

The MAC with conferencing encoders was introduced by 
Willems in [IJ, who also derived the capacity region for 
the discrete memoryless setting. Here we derive the capacity 
region for the Gaussian setting under average power con- 
straints. The achievability part is very similar to the one in 
lH). The converse, however, requires a novel tool first derived 
in |4| for proving that Gaussian distributions maximize certain 
mutual information expressions under a Markovity-constraint. 
For such maximization problems the traditional approach of 
proving the optimality of Gaussian distributions by employing 
the Max-Entropy Theorem |3j Theorem 12.1.1.] or a condi- 
tional version thereof E] fails. The reason is that replacing a 
non-Gaussian vector satisfying the Markovity condition by a 
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Gaussian vector of the same covariance matrix may result in 
a Gaussian vector that violates the Markovity condition. 

We also consider an additional scenario where the received 
signal is also corrupted by an additive Gaussian interference 
sequence which is non-causally known to both encoders but 
not to the decoder We show that, even though the decoder 
is non-informed, the interference can be perfectly canceled 
in the sense that the capacity region of the setting without 
interference is achievable also in this setting with interference. 

We next describe the channel models more precisely and 
then proceed to state our results. 

The goal of the transmission is that Transmitters 1 and 
2 convey their messages Mi and M2 to the receiver The 
messages Mi and M2 are assumed to be independent and 
uniformly distributed over the sets Mi = {1, . . . , [e"^ij } 
and A^2 = {1: ■ • ■ j [e"^^J}. Here n denotes the block-length, 
and Ri and R2 denote the rates of transmission in nats per 
channel use. 

Prior to each block of n channel uses, the two encoders 
hold a conference, i.e., they exchange information over k uses 
of the pipes. The pipes are assumed to be 

• perfect in the sense that any input symbol to a pipe is 
available immediately and error-free at the output of the 
pipe; and 

• of limited throughputs C12 and C21, in the sense that 
when the k inputs to the pipe from Transmitter 1 to Trans- 
mitter 2 take values in the sets Vi.i, . . . , Vi^k and the k 
inputs to the pipe from Transmitter 2 to Transmitter 1 
take values in the sets V2,i, . . . , V2,a;, then 

k k 

^log|Vi,,| < nCi2 and ^log|V2,f| < nC2i. (1) 

1=1 1=1 
Here and throughout all logarithms are natural logarithms. 

Note that the communication over the pipes is assumed to be 
held in a conferencing way, so that the £-th inputs Vu € Vi^i 
and V21 G V2.1 can depend on the respective messages as well 
as on the past observed pipe-outputs: 

Vi.t = fij{Mi,V2.i,...,V2x-i), (2) 
V2.i = hl{M2:Vi.i,...,Vi^-i), (3) 

for some given sequences of encoding functions {fi,i]^^i and 

{f2.i}\=i where 

Ji.t : Ml X V2.1 X ... X V2,£-i — > Vi,^ (4) 
J2.t ■■ M2 X Vi.i X ... X Vix-i — > V2,e- (5) 



We define an {n,Ci2,C2i)-conference to be the collec- 
tion of an integer number k, two sets of input alphabets 
{Vi,i, . . . , Vi^fc} and {V2,i, . . . , V2.k}, ™d two sets of encod- 
ing functions {/14, . . . , fi^k} and {/2,i, . . . , f2,k} as in dUi 
and dU, where n, C12, C21, A:, and the sets {Vi^i, . . . , Vi^fc} 
and {V2,i, . . • , V2,fc} satisfy ([B- 

After the conference, Transmitter 1 is cognizant of the se- 
quence V2 = (V2.1, . . . , V2,fc) and Transmitter 2 is cognizant 
of the sequence Vi = (Vi.i, . . . , Vij.). The channel input se- 
quences Xi = {Xi^i, Xi_„) and X2 = (^2,1, . . . , X2,„) 



can then be described with encoding functions ip 
as 

(«) 

- , J") 



(«) 



and 



(«) 



X^^^\">{Ah,Y2), 



^r(A^2,Vi), 



where 



(«) 

(n) 
V2 



Ml X V2,i X . , 
M2 X Vi 1 X . , 



X V2^k 

X Vi.fc 



(6) 
(7) 

(8) 
(9) 



Additionally, we impose an average block power constraint 
on both channel input sequences: 



n 



,t=i 



(10) 



The multiple-access channel is described as follows. For 
given discrete-time t and channel inputs xi^t,X2,t S the 
time t channel output Yt is 



Yt 



Xl,t + X2A + Zf, 



(11) 



where {Zt} models the noise corrupting the channel and is 
given by a sequence of independent and identically distributed 
(IID) zero-mean Gaussian random variables of variance cr^ > 
0. 

Based on the output sequence Y = (Yi,...,y„) the 
decoder applies a decoding function 0^"^ 

M"^A^ixA^2, 
to produce the message estimates Mi and M2, i.e.. 



(12) 



(13) 



An error occurs whenever (Mi, M2) ^ {Mi, M2). 

A rate pair (_Ri,i?2) is said to be achievable over the 
Gaussian MAC with conferencing encoders if there exist a 
sequence of {(n, C12, C2i)}-conferences, two sequences of 
encoding functions {Vi^\v2^''} as in (H) and (|9]l, and a 
sequence of decoding functions {t/)*^"^} as in (fT2] l such that 
the probability of error tends to as the block-length n tends 
to infinity, i.e.. 



lim Pr 



(Mi,M2) ^ {Mi,M2) 



= 0. 



(14) 



The capacity region C is defined as the closure of the set of 
all achievable rate pairs. 

We also consider an extension of the setting at hand in 
the sense of Costa's "Writing on Dirty Paper" channel |2|. 



Thus, we assume an additional additive Gaussian interference 
sequence which is non-causally known to both transmitters 
but not to the receiver There exist two different scenarios 
one could envision. A scenario where the transmitters learn 
the interference sequence before the conference, and thus the 
inputs to the bit-pipes can depend also on the interference; 
or a scenario where the transmitters learn the interference 
only after the conference. It turns out that the presented 
results do not depend on which of the two scenarios is 
considered. The converse we present holds also for the setting 
where the transmitters know the interference already before 
the conference; and the encoding scheme applies also to the 
setting where the transmitters know the interference only after 
the conference. In the following we will focus on the setting 
where the transmitters learn the interference sequences after 
the conference. 

For the setting with interference we need to modify the 
definitions of the channel in ( fTTT l and the encoding functions 
in (O and (|9]i; the decoding function remains as in (fT2l i. 

For given inputs xi t and X2 t the channel output is given 

by 



Yt = xit+X2,t + St + Zt, 



(15) 



where the noise sequence {Zi} is defined as before, and where 
the interference sequence {St] is an IID zero-mean Gaussian 
sequence of variance Q and independent of the noise sequence 
and of the messages. 

Denoting the interference sequence by S = . . . ,Sn), 
the channel input sequences Xi and X2 are described as 

Xi =^("i(Mi,V2,S), 
X2 = ^^",^(M2,Vi,S), 

for some encoding functions <Pi"p, ^2% of the form 



(n) 

Vi,iF : 

(") 
<^2,iF : 



Ml X V2A X 
M2 X Vi.i X 



X V2,fc X 

X Vi.fc X 



The input sequences are subject to the power constraints ( fTOl l. 

The probability of error, achievable rate pairs, and capacity 
region Cif for this new setting are defined as before. 

II. Main Results 
Definition 1: Define the region 



Cc 



U 1(^1,^2): 

0</3i,/32<l I 

Ri < ilog(^i + ^^ +C12, 

i?2<^l0g(^l + ^) +C21, 

-Ri + -R2 < :r log 1 + 2 I + C^i2 + C21, 



p^p^l] L , + ^2 + 2v/PiP2/3l/32 \ 1 

i?l + i?2 < 2 log IH I >, (16) 



where Pi — 1 — (3i and = 1 — /32- 

Theorem 1: The capacity region C of the Gaussian MAC 
with conferencing encoders is equal to Cg, 



C=Cr 



(17) 



Remark 1: The main step in the proof (see Section BlI-Al i 
is to show that under the Markov condition Xi —o—U -~o—X2 
the region 7?.conf(^i, C^, -^^2) (Definition |2]) is maximized by 
choosing jointly Gaussian distributions on X2). A 

subset of the mutual information expressions which char- 
acterize the region 7?.conf(-''^ij C^j -''^2) can be found in the 
characterization of an achievable region for the MAC with 
perfect feedback proposed by Cover and Leung |6| and in 
the characterization of the capacity region of the MAC with 
common messages derived by Slepian and Wolf |7|. Again, in 
both characterizations the triple Xi^>—U^>—X2 is required 
to be Markov, and the same tools as in Section UlI-AI can be 
used to prove that for Gaussian channels also the Cover-Leung 
region and the Slepian-Wolf region are maximized by choosing 
jointly Gaussian distributions on Xi^>—U^>—X2 iH). 

Theorem 2: The capacity region Cif of the Gaussian MAC 
with conferencing encoders and additive Gaussian interference 
sequence non-causally known at both encoders equals the 
capacity region C of the setting without interference 



CiF — C . 



(18) 



Note that Costa's result Q on "Writing on Dirty Paper" and 
Gel'fand and Pinsker's result fSl on "Multi-access Writing on 
Dirty Paper" are special cases of Theorem |2] 

III. Proof of Theorem[T] 
The achievability of Cg, i.e., 

CgQC, (19) 

follows by applying the scheme described in |[T1 with a 
Gaussian input distribution. The details are omitted. 

A. Converse 

To prove the converse, i.e.. 



(20) 



we first outer bound C by Cout (Lemma[T]). The converse is then 
established by showing Com = Cg. To this end, in Lemma|2]we 
express the region Cg in a similar form to Cout, i-C-, as a union 
of regions where the union is taken over certain distributions 
satisfying a Markov condition and power constraints. We then 
notice that Cout and Cg differ only with respect to the set 
of distributions over which the unions are taken (see (ISTT i 
and (l23Tl): for Cout the union is taken over all distributions 
satisfying the Markov condition and the power constraints, 
and for Cg the union is only over those that are Gaussian. We 
conclude the proof by showing that for Cout it is sufficient to 
take the union only over Gaussian distributions (Lemma [3]). 



Definition 2: For a given distribution pXlUX2{^^ 'i on the 
random triple Xi, [/, X2, define 



7?.Conf(-''^l, U, X2 

= [{Ri,R2) ■■ 



Ri <IiXi;Y\X2U) + Ci2, 

R2 <I{X2;Y\XiU) + C2i, 

R1+R2 <I{XiX2:Y\U) + Ci2 + C2u 

R1+R2 <I{XiX2;Y) }, 

where the mutual informations are computed 
with respect to the law PuXiX2Y{u,xi,X2,y) = 
PXtUX2ixi,u,X2)p{y\xi,X2), where p{y\xi,X2) denotes the 
channel law. 

Definition 3: Define the region 



Cout — [J Ti-Coni{Xi,U, X2) 

Xi^-U-^-X-i 
E[xl]<Pt, E[XI]<P2 



(21) 



where the union is over all joint distributions (not necessarily 
Gaussian) for which Xi^>— [/^>— ^2 is Markov and for which 
E [X'l] < Pi and E [X'i] < P2 

Lemma 1: The region Cout is an outer bound on the capacity 
region of the Gaussian MAC with conferencing encoders. 



C C Cout- 



(22) 



Proof: Requires only a slight modification of the outer 
bound in |1J to account for the power constraints. ■ 
Lemma 2: The region Cg in ( fT6l ) can be expressed as 

Cg= U nconf{X^,U^,X^), (23) 

xf^-u^^-x§ 
e\{x?Y]<p,.e\{x§Y]<P2 



where the superscript Q is used to indicate that the union is 
taken only over Gaussian Markov distributions satisfying the 
second moment constraints. 

Proof: Follows by evaluating the various mutual informa- 
tion terms in the definition of 7?.conf for Gaussian distributions 
on Xi^^U^^X2. m 
The right hand-sides of (ISTT l and (l23T l differ only with 
respect to the set of distributions over which the unions are 
taken. Therefore, in order to conclude the proof of the converse 
(I20I 1. by Lemma [U and Equations (l2Tl l and (|23] |. it suffices 
to show that there is no loss in optimality if the union in 
(ISTT i is taken only over Gaussian Markov triples fulfilling 
the second moment constraints. This is established by the 
following Lemma |3] 

Lemma 3: For any Markov triple Xi~o—U —0—X2 fulfilling 
E [Xf] < Pi and E < P2, there exists a Gaussian 

Markov triple X^*-o-V'^*-o-X^* fulfilUng the power con- 
straints E^fXpi^l < Pi and E[(xf*)^j < P2, such that 



nconf{Xi,U, X2) C ncondX^*,V^*,X^*). (24) 

We postpone the proof of Lemma [3] and first state a sequence 
of definitions and lemmas. 



Lemma 4: For any (not necessarily Markov) random triple 
X\ , U, X2 of finite second moments, 

7^Corf(^i, U, X2) C 7^conf(^f , U^, Xi) 

where {Xf,U^,X^) is a centered Gaussian vector whose 
CO variance matrix is equal to that of (Xi, J7, X2). 

Proof: Follows by a conditional version of the Max- 
Entropy Theorem [3, Theorem 12.1.1.], see also fSl. ■ 
Definition 4: Define JCg as the set of 3 x 3 positive semi- 
definite matrices 



Proof of Lemma \3\ Let 



K 



fell fei2 fcl3 
fcl2 ^22 k-23 
fcl3 ^23 ^33 



(25) 



satisfying one of the two conditions 

1) /c22 and fci3fc22 = ^12^23; 

2) /C22 = A:i2 = /ci3 = fc23 = 0. 

Lemma 5: A Gaussian triple is Markov if, and only if, its 
covariance matrix is in K.g. 

Proof: The "if" direction follows because the law of 
a Gaussian triple is fully characterized by its mean and its 
covariance matrix and by noting that for any covariance matrix 
K e JCg we can construct a Gaussian Markov triple of 
covariance matrix K. The proof of the later is omitted. 

For the proof of the "only if" direction we assume a 
triple A^—B^>—C which forms a Markov chain in this 
order. We distinguish between two cases: Var(i3) = and 
Var(i3) 7^ 0. If Var(i3) = 0, then B is deterministic 
and the Markov chain A^>—B^>—C implies that A and C 
are independent. The covariance matrix of A,B,C is then 
diagonal, and Condition |2]l is satisfied. If Var(i3) ^ 0, then 
define Aq ^ A- E[A] , Bq = B - E[B], and Co = C-E[C] 
and compute 

E[AoCo] - E[E[^Co|So]] = E[E[Ao|So] E[Co|Bo]] 

E[AoBo]E[BoCo] 



E[AoBo] EjBoCo]^- 



(26) 



Var(i?o) "Var(i?o) 
Here, the second equality follows by the Markovity and the 
third equality by the Gaussianity. By multiplying (|26] | with 
E[Bo] = Var(i3) Condition [B is obtained. 

Lemma 6: Consider a Markov triple Xi- 
Xi and X2 of finite second moments. Let 



-U^>-X2 with 



V = E[Xi\U] - E[Xi] . (27) 

Then, the covariance matrix of the triple {Xi, V, X2) is in ICg, 
and 

7^co„f(^l, U, X2) C 7^co^f(Xl, 1/, X2). (28) 

Proof: The inclusion (l28T l follows by the following two 
observations. Exchanging [/ by a deterministic function of U 
increases all mutual information expressions in 7?.conf which 
are conditional on U. And changing U does not change 
the joint distribution of Xi, X2, and hence the unconditional 
mutual information expression remains the same. 

That the covariance matrix of the triple {Xi,V,X2) is 
in JCg follows because (|27| i and the Markov condition 
Xi^y-U^:>-X2 imply that 

Cov[V, X2] = Cov[Xi, X2] and Cov[V, Xi] = Var{V) . ■ 



V^E[X,\U]-E[X,], 



(29) 



X2* to be zero-mean jointly 



and define the triple Xf*,V^*,X^* 
Gaussian with the same covariance matrix as the triple 
Xi ,V,X2- To conclude the proof we shall show that 
X^* ^y-V^* ^>-X^* forms a Markov chain and that Con- 
dition dH is satisfied. That the triple Xf*-<^V^*^X^* 
forms a Markov chain follows by the Gaussianity of 
Xf* ,V^* , X^* and by Lemma |5l because by Lemma 
the covariance matrix of Xi , V, X2 is in JCg and thus, by 
construction, also the covariance matrix of X^* ,V^* , X2* 
is in JCg. Note that the triple Xi,V,X2, even though its 
covariance matrix is in ICg, does not necessarily form a 
Markov chain in this order, because it is not restricted to be 
Gaussian. 

That Condition (l24l i is satisfied follows by the following 
sequence of inclusions: 

7^Conf(^l, U, X2) C 7^Conf(^l, V, X2) 

c7^conf(^f^l^^^^2^*), (30) 

where the first inclusion follows by Lemma |6] and the second 
inclusion follows by Lemma |4] ■ 



IV. Proof of Theorem[2] 



The "converse" 



CiF CC 



(31) 



follows because C outer bounds the capacity region of the 
channel with interference even when the interference is also 
known at the receiver It remains to prove the "direct part" 



CQCu 



(32) 



i.e., that every rate pair in C is achievable in the presence of 
interference. This follows by Lemmas [T] and [8] ahead. 
Definition 5: Define the region 



7^Ach= IJ \{Ri,R2): 

0</3i:/32<l I 



^1 < log 1 



1 



f3iPi 



(33) 



2 \ PiPi + fj2P2 + (T-^ / 



i?2<^log('l + ^) +C21, 



i?2 < ^ log ( 1 



I32P2 



1, {^/KPi + VNP2? 

2 /3lPl+/32P2 + ^2 



(35) 



(36) 



i?i+i?2 < ^logfl + ^i^^-^2^ ) +C12 + C21, (37) 



i?i + i?2 < 2 log H I } ' (38) 

where f3i — 1 — Pi and ~ 1 — ^2- 

Lemma 7: The capacity region Cw includes 7?.Ach: 

CiF 2 7^Ach. (39) 

Proof: See Section HV^ ■ 
Lemma 8: The achievable region 7?.Ach equals the capacity 
region C of the Gaussian MAC with conferencing encoders, 

A. Coding Technique Achieving TZacH 

In this section we sketch a coding technique that achieves 
the region TZaoi- The analysis is omitted. 

The two transmitters first create a common message by 
communicating over the pipes as in |1|. Thus, after the 
conference Transmitter 1 is cognizant of the common message 
and of an independent private message. It allocates power 
(1 — f3i)Pi to the common message and power PiPi to the 
private message. Similarly for Transmitter 2. 

The coding technique involves time-sharing between two 
schemes. Both schemes apply successive decoding at the 
receiver, where the receiver first decodes the common message 
followed by the private messages. But they differ in the de- 
coding order of the private messages. We describe the scheme 
where the decoding of the common message is followed 
by the decoding of the private message of Transmitter 1 
and only thereafter by the decoding of the private message 
of Transmitter 2. The other scheme where the decoding of 
the common message is followed by the private message of 
Transmitter 2 is analogous. 

We first describe the encoding of the common mes- 
sage. Before transmission begins, the transmitters agree 
on a (single-user) dirty-paper code for power Po — 

(^(1 - l3i)Pi + 1/(1 - ^2)^*2) , noise variance (/3iPi + 
/32P2 + f^), and interference S. Transmitter 1 encodes the 
common message using this dirty -paper code and scales the 

resulting sequence by ^^^^ . Transmitter 2 encodes the 
common message with the same code, but scales the resulting 

sequence by (The channel will coherently com- 

bine the two sequences.) 

Independently of the common message, the transmit- 
ters encode the private messages. Transmitter 1 encodes 
its private message using a dirty-paper code for power 
(3iPi, noise variance + /32-P2, and interference Si = 

( 1 iV^Ef^^^^l_] s. Transmitter 2 en- 

codes its private message with a dirty-paper code for 
power (32P2, noise variance a^, and interference S2 ^ 
^1 — ff^p^+/;fp3+(T^ ^ Si. Each transmitter sends the sum of 
the two sequences produced for the common message and for 
its private message. 



The receiver performs successive decoding. It starts by 
decoding the common message based on nearest neighbor 
decoding as in |9], while treating the sequences which the 
transmitters produced for the private messages as additional 
noise. Then, the receiver subtracts (or "strips off") the decoded 
common-message codeword from the channel outputs and pro- 
ceeds to decode the private message of Transmitter 1 . (Here, 
the "common-message codeword" is not the resulting sequence 
of the dirty-paper code, but the codeword in the bin of the 
common message which was selected during the encoding 
procedure of the dirty-paper code.) To decode the private 
message of Transmitter 1, the receiver again uses nearest 
neighbor decoding and treats the sequence which Transmitter 2 
produced for its private message as additional noise. Finally, it 
subtracts the decoded Transmitter 1 -private-message codeword 
and decodes the private message of Transmitter 2. 

Remark 2: To encode the different messages the transmit- 
ters use dirty-paper codes for scaled versions of the interfer- 
ence S. The reason for this is that the codewords that the 
receiver subtracts depend on the interference sequence S, and 
the resulting channel seen in subsequent decoding phases is 
interfered by a scaled version of S. 

Remark 3: Our coding scheme is different from Willems's 
scheme |[l] (for the setting without interference). In his scheme 
the transmitters apply superposition coding and the receiver 
applies joint decoding. Our approach has two advantages: It 
simplifies the analysis, and it achieves the same result also if 
the noise sequence {Zt} is not IID Gaussian but any arbitrary 
ergodic process of second moment a^. However, our approach 
leads to the additional constraints ( l34l l and ( |36] |. Fortunately, 
by Lemma [8] these additional constraints don't shrink the 
resulting region. 

Acknowledgment 

We would like to thank Stephan Tinguely, i. Emre Telatar, 
and Tsachy Weissman for helpful discussions. 

References 

[1] F. M. J. Willems, "The discrete memoryless multiple access channel 
with partially cooperating encoders," IEEE Trans, on Inform. Theory, 
vol. 29, no. 3, pp. 441^45, Nov. 1983. 

[2] M. H. M. Costa, "Writing on dirty paper," IEEE Trans, on Inform. 
Theory, vol. 29, no. 3, pp. 441^45, Nov. 1983. 

[3] T. M. Cover and J. A. Thomas, "Elements of information theory", 2nd 
edition. Wiley publication. 

[4] V. Venkatesan, "Optimality of Gaussian inputs for a multi-access achiev- 
able rate region". Semester Thesis, ETH Zurich, Switzerland, June 2007. 

[5] J. Thomas, "Feedback can at most double Gaussian multiple access 
channel capacity," IEEE Trans, on Inform. Theory, vol. 33 no. 5, Sep. 
1987. 

[6] T.M. Cover and C.S.K. Leung " An achievable rate region for the 

multiple-access channel with feedback," IEEE Trans, on Inform. Theory, 

vol. 27 no. 3, pp. 292-298, May. 1981. 
[7] D. Slepian and J.K. Wolf, "A coding theorem for multiple access 

channels with correlated sources," Bell System Tech. J. vol. 52, pp. 1037- 

1076, Sept. 1973. 

[8] S. Gel'fand and M. Pinsker, "On Gaussian channels with random 
parameters," in Proc. ISIT 1984, Taskent, U.S.S.R., Sep. 1984, pp. 247- 
250. 

[9] A.S. Cohen and A. Lapidoth, "The Gaussian Watermarking Game", 
IEEE Trans, on Inform. Theory, vol. 48, no. 6, June 2002. 



